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ABSTRACT 

We show that a large-scale, weak magnetic field threading a turbulent accretion disk tends to be 
advected inward, contrary to previous suggestions that it will be stopped by outward diffusion. The 
efficient inward transport is a consequence of the diffuse, magnetically-dominated surface layers of 
the disk, where the turbulence is suppressed and the conductivity is very high. This structure arises 
naturally in three-dimensional simulations of magnetorotationally unstable disks, and we demonstrate 
here that it can easily support inward advection and compression of a weak field. The advected field 
is anchored in the surface layer but penetrates the main body of the disk, where it can generate strong 
turbulence and produce values of a (i.e., the turbulent stress) large enough to match observational 
constraints; typical values of the vertical magnetic field merely need to reach a few percent of equipar- 
tition for this to occur. Overall, these results have important implications for models of jet formation 
which require strong, large-scale magnetic fields to exist over a region of the inner accretion disk. 

Subject headings: accretion, accretion disks — galaxies: jets — magnetic fields — MHD — X-rays: 
binaries 



1. INTRODUCTION 

Early theoretical work on accretion disks argued 
that a large-scale magnetic field (of, for exam- 
ple, the interstellar medium) would be dragged in- 
ward and greatly compressed by the a ccreti n g plasma 
(iBisnovatyi-Kogan fc RuzmaikinI Il974l . 119761 : iLovelacd 
Il976l ) . Figure [1] illustrates this concept by showing a 
sketch of an ordered magnetic field threading an accre- 
tion disk, in which inward advection has caused the mag- 
netic field lines to bunch together into an "hourglass" 
shape. This was thought to be a simple mechanism for 
generating dynamically significant fields in the inner disk. 

In the present paper, we revisit this issue, build- 
ing off the recent work of iBisnovatvi-Kogan &: Lovelacd 
12007). Our motivation for doing so is that in the in- 
tervening years, the early theoretical arguments have 
been challenged. More detailed models of turbu- 
lent disks suggested that a large-scale, weak mag- 
netic field such as that show n in Figure [T] in fact 
will diffuse outward rapidly (Ivan Ballegooije 3 Il989l : 
iLubow. Papaloizou. fc Pringld ll994J l if the turbulent 
magnetic diffusivity and turbulent viscosity are of sim- 
ilar order of magnitude, as they are expected to 
be (lParkeilll97H iBisnovatvi-Kogan fc RuzmaikinI Il976l : 
ICanuto fc Battaglial Il988l) — the turbulence responsible 
for driving the accretion also leads to enhanced recon- 
nection of the large-scale radial field across the thickness 
of the disk, thereby causing the vertical field to diffuse 
away. This cast doubt on the idea that weak fields could 
be dragged inward and compressed by advection. At 
the same time, it was known that the angular momen- 
tum loss to magnetohydrodynamic (MHD) outflows from 
a disk threaded by a sufficiently strong large-scale field 
could more than offset the outward diffusion and lead to 
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a rapid, implosive i ncrease of the field in the c entral re- 
gion of the disk (Lovela ce. Romanova. fc Newm an 1994). 
However, it seemed to be the case that growth of a strong 
magnetic field "from scratch," due to continual advection 
of a weak field, was impossible in a thin disk. Although 
this conclusion has been occasionally challenged (e.g., 
lOeilvie fc Lhdoll200l[l . it is still generally accepted, which 
has led to the recent suggestion that special conditions 
(extremely nonaxisymmetric regions of strong field in an 
otherwise weakly-magnetize d disk) are required for t he 
field to be advected inward llSpruit fc Uzdenskvl l2005^. 

At the same time, recent three-dimensional MHD 
simulations have been performed that allow this issue to 
be addressed computationally. These simulations resolve 
the largest scales of magnetorotational turbulence and 
therefore self-consistently include the turbulent viscosity 
and diffusivity (without having to prescribe their values 
a priori). Most simulations performed to date have 
investigated conditions in which the accreting matter 
does not contain any net magnetic fiux and where 
no magnetic field is supplied at the boundary of the 
computational domain. However, in one simulation, 
weak poloidal flux injected at the outer boundary was 
clearly observed to be dragged into the central region 
of the disk, lea ding to the buildup of a strong central 
magn etic field (jlgumenshchev. Naravan. fc Abramowicd 
I2003D . A similar process, albeit transient, may occur in 
simulations without a net magnetic flux; there, radial 
stretching of locally poloidal field lines in the initial 
configuration often leads to large-scale p oloidal fields 
and i e t structures in the inne r disk (e.g.. iHirose et al.l 
I200I iDe Villiers etrall 1200^ iHawlev fc KrolikI l2n(ift 
see also the discussio n in llgumenshchev et al.l l2003l 



and iMcKinnev fc NaravanI 20071 and espe cially 



the simulations of iMcKinnev fc Gamnii^ l2004l and 
iBeckwith. Hawlev. fc Krolild I2007L which explore the 
effect of different initial field geometries on the forma- 
tion of jets). The extent to which any of the advection 
of magnetic field lines seen in numerical simulations 
requires the presence of a thick disk or nonaxisymmetric 
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Fig. 1. — Sketch of the magnetic field threading an accretion 
disk, showing the increase of the field due to its assumed inward 
advection with the gas (as proposed in early theoretical models). 



conditions is unclear. 

In light of these numerical results, we return to the 
question of inward advection of magnetic fields in this 
paper, allowing for the possibility that the disk is thin 
and axisymmetric and asking once again whether advec- 
tion of a weak field is possible under these conditions. 
The mechanisms we discuss here can occur in sufficiently 
ionized regions of any accretion disk; although they are 
perhaps most widely applicable to disks around black 
holes (where the large-scale magnetic field arises entirely 
within the accreting plasma), they are relevant for disks 
around many other types of accreting objects as well. 

The organization of this paper is as follows. In ^ we 
analyze the advection of a large-scale field in an accre- 
tion disk and point out the importance of the vertical 
structure of the disk, which was not taken into account 
in most p revious studies. Based on a n earlier suggestion 
(|Bisnovatvi-Kogan fc Lovelacell2007D . we show that the 
thin, highly conducting surface layer of the disk, where 
turbulence is suppressed, allows a large-scale magnetic 
field to be advected inward and compressed. In Sj3l 
argue that the resulting magnetic flux through the main 
body of the disk (due to the large-scale field being ad- 
vected inward) can produce values of the turbulent a 
parameter that are in accord with observational data. 
This is in contrast with numerical simulations of turbu- 
lent disks without a net imposed magnetic flux, which 
are unable to generate large enough turbulent stress. Fi- 
nally, in 21 we derive detailed conditions on the fleld 
strength, geometry and ionization fraction that are re- 
quired for the field to be advected inward and show that 
these are typically weak constraints. Conclusions of this 
work are summarized in fj5l 

2. MAGNETIC FIELD ADVECTION AT THE SURFACE OF 
AN ACCRETION DISK 

The evolution of the magnetic field B in an accretion 
disk (averaged over the short timescales of the turbu- 
lence) is assumed to be described by the induction equa- 
tion, 

— = V X (v X B - 7/V X B) , (1) 

where v is the plasma velocity, 77 — c^/(47rcr) is the mag- 
netic diffusivity, c is the speed of light, and a is the con- 




FlG. 2. — Sketch of the disk and instantaneous poloidal magnetic 
field lines considered in this work. The toroidal field component 
is not shown. The inset shows a rough illustration of the vertical 
profile of the conductivity a (z) in units of the coronal value a (h). 
At the base of the corona (z = h), turbulence is suppressed and the 
conductivity is very high; therefore, if the material in this region 
advects inward with the main body of the disk, the large-scale 
magnetic field will be advected inward as well. 

ductivity.^ 

We assume a disk with half-thickness _ff < r in cylin- 
drical coordinates. The main body of the disk is turbu- 
lent, and we take the effective diffusivity to be ry ~ v, 
where ly is the turbulent viscosity. The turbulence is 
widely th ought to be due to the magnetorotational in- 
stability iB albus fc Hawle^ [l99ll . llQQSt IVelikhovl HOSOI : 
I Chandraie khar I960), which roughly occurs when the 
magnetic energy density is less than the thermal energy 
density. We therefore assume a weak magnetic field such 
that this condition holds in the main body of the disk. 

However, the time-averaged magnetic field is not ex- 
pected to vary strongly across the disk thickness, owing 
to the buoyancy of the field and the condition V • B = 
(in more physical language, the field is not influenced 
by the vertical gravity that keeps disk material con- 
fined near the equatorial plane). Thus, the mass den- 
sity of the gas will typically decrease with height z more 
rapidly than the time-averaged magnetic field strength, 
and at a height ~ H above the midplane, the magnetic 
energy density will become strong enough compared to 
the thermal energy density that turbulence will be sup- 
pressed. The boundary between the turbulent and non- 
turbulent regions is likely to be "fuzzy" owing to the leak- 
age of some magnetic flux throu gh the disk surface (e.g., 
iGaleev. Rosner. fc Vaianalll979f ). but at a certain height, 
the plasma will become completely nonturbulent. In this 
paper, we will use the terms "base of the nonturbulent 
region" and "surface layer of the disk" interchangeably; 
however, it should be noted that we are explicitly defln- 
ing these regions to be above the boundary layer and 
therefore fully a part of the nonturbulent corona (see 
Figure [21). 

* Note that in a turbulent disk, there can also be an additional 
term in equation JTJ that we have not included here, which would 
represent the contribution of a turbulent dynamo to the growth 
of the large-scale field. We ignore this term because we are only 
interested in the growth of magnetic field due to advection, and 
therefore any local, dynamo-generated field that may be produced 
is "extra" to that which we discuss in this section. 
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This suppression of turbulence above a weakly- 
magnetized disk has been o bserved in a vari- 
ety of MH D simulati ons (e g.. iMiller fc Ston 
De Villiers. Hawlcy, & Krolild l2003t iHirose et~all 



2000 



2004 



McKinnev fc Gammid |2004|; iFromang fc Nelson' 2006*) 



inclu ding those with radiation |Hirose. Krolik. fc Stone 
[200l) and even those in which t he radiation pressure is 
comp arable to the gas pressure (jKrolik. Hirose. fc Blaed 
120071) . Howeve r, MHD simu lations of fully radiation- 
dominant disks (|TurneHl2004j) are less clear, and the ap- 
plicability of our work in this case requires further anal- 
ysis. Nonetheless, even above a radiation-dominated re- 
gion of the disk, we expect that the turbulence will be 
suppressed in many situations; we discuss this issue fur- 
ther in 21 

Figure [2] shows a schematic drawing of the considered 
geometry. The lack of turbulence at a height h ^ H 
near the disk surface causes this layer to become highly 
conductive; the diffusivity will decrease from its turbu- 
lent value in the main body of the disk (77 ~ ~ 10^^ 
cm^ s~^ for typical parameters) to the Spitzer value as- 
sociated with electrons scattering off of ions, given by 

200 (Ts/ keV)"^/^ cm^ s-\ where T, is the surface 



^75 



temperature.^ This suggests that the second term on the 
right hand side of equation ([1]) can be ignored in the up- 
per disk layers. Specifically, the relative importance of 
the two terms (advection compared to diffusion) at any 
point in the disk is determined by the local magnetic 
Reynolds number Rem = Hur/t]^ where Ur is the local 
radial speed and H is the relevant length scale (here we 
make the reasonable assumption that the time-averaged 
magnetic field does not vary significantly in the radial di- 
rection on l ength scales shorter than H ). We can there- 
fore use the IShakura fc SunvaevI ()1973l ) disk solution to 
find that a typical value of the magnetic Reynolds num- 
ber at the surface of the disk, where turbulence is sup- 
pressed, is given by 

Rem - 10* «m5/8m3/8f-5/«/*"'^' {lO^H/rfUs. (2) 
Here, a < 1 is the dimensionless "viscosity" parameter 
in the main body of the disk (i.e., the magnitude of the 
turbulent stress divided by the thermal pressure), m is 
the mass of the central object in solar masses, m is the 
accretion rate in units of the Eddington luminosity di- 
vided by the speed of light squared, f is the radius in 
units of the Schwarzschild radius, /, < 1 is a dimen- 
sionless function of f that depend s on the stress at th e 
inner boundary of the disk (e.g.. lAgol fc KrolikI [200I . 
and Us is the ratio of the radial speed at the disk surface 
to that in the main, turbulent body of the disk. This last 
term can be smaller than unity, but not small enough to 
prevent the conclusion that, typically. Rem ^ 1 at the 
surface of the disk, and diffusion of the magnetic field 
can be neglected. This is in contrast with the main body 
of the disk, where Rem ~ H/r (assuming the magnetic 
field is not strong enough to affect the accretion speed) 
and diffusion o f the magnetic field therefore dominates 



over advection ( | Lubow et al 



iHevvaerts. Priest, fc Bardoul 



[T99I iLovelace et al.lll99l 



T9M ). 



^ In actuality, other non-ideal MHD effects, in particular Ifall 
electromotive forces, m ay be more important than the Spitzer dif- 
fusivity, as we show in i|4.3l However, the Hall effect usually does 
not oppose inward advection of the magnetic field, and even if it 
does, it will not be important unless the field is very weak. 



We can easily demonstrate that advection in the sur- 
face layer of the disk is able to support the overall growth 
of magnetic field. If we integrate equation (jT]) over a cir- 
cular surface r < tq that covers the top side of the disk 
(z = /i, where h is the height at which turbulence is first 
suppressed), we can take 77 « 0, and Stokes' theorem 
therefore implies that 

=rQ dcj} {VzhBrh - VrhBzh) \r=ro , (3) 

where $p is the poloidal magnetic flux through this sur- 
face and the h subscript indicates that the quantity is 
evaluated aX z = h. If the right hand side of this equa- 
tion has the same sign as $p, the magnetic flux interior 
to radius tq will grow. 

Assuming axisymmetry (or, alternatively, treating sub- 
sequent quantities as being appropriately averaged over 
azimuth) and taking z > for definiteness, the condition 
for magnetic fiux growth in equation ([3]) simplifies to 



(4) 



provided that magnetic field with the appropriate polar- 
ity is available to accrete. 

Although turbulent stress cannot contribute directly to 
the accretion at z = h, coupling between the main, tur- 
bulent body of the disk and the surface (as well as angu- 
lar momentum loss to a wind or jet) will tend to produce 
Vrh < 0. Also, simulations indicate that for an MHD 
outfiow or jet (vzh > 0), the magnetic field structure has 
Brh/B^h > (|Ustvugova et al.l ll999, 2000). Thus, equa- 
tion ^ is in general likely to be satisfied, a point which 
we will discuss more rigorously in §4] 

In summary, our arguments in this section are a sim- 
ple consequence of the fact that magnetic fields are sus- 
tained by the fiow of current, not the fiow of mass. In 
order to prevent inward advection of magnetic fields, tur- 
bulent diffusion must oppose advection throughout the 
entire inward-accreting portion of the disk, so that no 
currents are allowed to accrete inward. Even a small 
sliver of nonturbulent (i.e., highly conducting) material 
that advects inward at the surface layer can support the 
magnetic field, even though it may only contain a small 
fr action of the disk's ma ss. A related issue was noticed 
bv lOgilvie fc Liviol ()2001l ), who argued (on mathematical 
grounds) that the relevant radial velocity for magnetic 
field advection is one that has been weighted by l/?7 and 
averaged over height. Here, we present a physical model 
for the behavior of ry with height and show that in a disk 
where the magnetic diffusivity is due to turbulence, the 
contrast between the diffusivity inside and outside the 
turbulent region is likely to be so sharp that the condi- 
tion for magnetic flux growth reduces to equation (U, 
which is satisfied in many parts of a typical accretion 
disk. 

3. EFFECT OF THE ADVECTED MAGNETIC FIELD ON 
THE TURBULENT a PARAMETER 

Given the apparent ease with which a large-scale mag- 
netic field can advect inward in an accretion disk, it is 
natural to consider the infiuence of this magnetic field on 
the disk dynamics. In particular, in this section we dis- 
cuss how the advected magnetic field might be expected 
to affect the turbulence in th e main body of the disk, em - 
bodied in the a parameter of lShakura fc SunvaevI (|1973f) . 
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Since the turbulence is thought to be magnetic in nature 
and in particular due to the magnetor otational instability 
(MRI; iBalbus fc Ha^levl [l99lL Tl99l . the efJect is likely 
t o be a significant one. 

iKing . Pringle. fc Livid ()2007| ) have recently pointed 
out that observationally-determined values of a, based 
primarily on studies of outbursts in dwarf novae and 
X-ray transients, tend to lie in the range ~ 0.1 — 0.4. 
Recalling that a is a measure of the turbulent magnetic 
stress scaled by the thermal pressure in the main body 
of the disk, it is clear that a significant amount of 
turbulent magnetic energy must exis t in these ac c retion 
disks during the outburst phase. iKing et all (|2007l ) 
noted a potential puzzle, which is that numerical 
simulations of the MRI in which the instability is 
allowed to develop entirely based on a local seed field 
(i.e., where there is no externally-imposed magnetic 
fiux through the computational region) tend to give 
saturation values of the stress that are much too 
small to match the observations, with typical values 
a ^ 0.01 regardless of the s treng th of the seed field 
(iHawlev. Gammie, fc Balbui 11996: iBalbus fc HawlevI 



19981 and note that iPessah. Chan, fc PsaltisI l2007l and 



Fromang fc Papaloizoul l2007l have shown that even 
these values may be significant overestimates, due to 
numerical resolution effects). 

However, in simulations with an externally-imposed 
vertical magnetic field, the turbulent stress due to 
the MRI depends critically on the seed field strength. 
In particular, a is found to increase with the net 
impos ed vertical field (iHawlev. Gammie. fc Balbud 
Il995t IBalbus fc HawlevI mM . Shearing box simu- 
lations suggest a rough empirical relationship a ~ 

27r {Pz,ext) , where Pz.ext is the ratio of the ther- 
mal pressure in the disk to the magnetic pressure of 
the externally-imposed vertical field (this is a sim- 
plified version of a more general equation found in 
Pessah et al.ll2007L which is based on MRI simulations by 



Sano et al.ll2004 but a lso agrees with the earlier results 
of IHawlev eranil995D . Thus, values of a ~ 0.1 - 0.4 
simply indicate a moderately strong (but still signif- 
icantly sub-equipartition) vertical field, perhaps with 
~ 250 - 4000. 

iKing et al] ()2007| ) ruled out this mechanism because 
they believed there was no obvious source for such an 
externally-irnposed field in real accretion disks (whereas 
iPessah et al.l l2007l believed it to be possible, but at- 
tributed the field to internally-generated MRI fluctua- 
tions, rather than an external source). However, if large- 
scale magnetic fields can be advected inward in a disk, 
and if, furthermore, these fields are maintained by cur- 
rents flowing in the nonturbulent surface layer of the 
disk, as we have argued, then MRI simulations with an 
externally-imposed vertical field are in fact the most rel- 
evant for comparing with real accretion disks. The net 
vertical field is anchored in the surface layer, and the 
main, turbulent body of the disk sees this field as a fixed, 
"externally-imposed" seed field. 

We therefore conclude that a typical value for the verti- 
cal magnetic field in an accretion disk undergoing an out- 
burst, based on observations, is ~ 2 — 6% of the equipar- 
tition field strength (i.e., Pz,ext is ^ 250 — 4000) and 
that there is no difficultly reconciling the observationally- 



determined values of a with numerical simulatio ns of 
the MRI. Also, as noted by iKing et all (|2007l ). the 
observationally-determined values of a are weighted av- 
erages over the entire accretion disk. Thus, it is certainly 
possible that advection could lead to much larger field 
strengths in a particular region. The field strengths in 
quiescent disks are similarly unconstrained. 

An interesting effect of the dependence of a on the 
strength of the large-scale magnetic field is that the vis- 
cous and thermal timescales in the disk (which depend 
inversely on a) should change with time, in response to 
the history of magnetic field advection. This process may 
explain some of the wide range of variability on many 
different timescales seen in accreting black holes, in par- 
t icular in the bright X-ray binary GRS 1915-1-105 (e.g., 
iBelloni et al . 2000) , where the various modes of variabil- 
ity seem to repeat in a semi-regular pattern over a period 
of months to years that has b een suggested to b e a sig- 
nature of magnetic processes ()Tagger et al.l [200^ . 

A large-scale magnetic field may be expected to have 
other effects on an accretion disk besides those discussed 
above. If the field is strong enough, it can begin to affect 
the disk dynamics directly (through removal of angular 
momentum via a wind or jet); we will consider the case 
where the advected field builds up to dynamically signif- 
icant values in a future paper. However, the important 
point we make in this section is that even when the large- 
scale field is dynamically weak, it can have a significant 
effect on the disk dynamics indirectly, through its influ- 
ence on the turbulent stress in the main body of the disk. 



4. DETAILED ANALYSIS OF THE CONDITIONS FOR 
MAGNETIC FIELD ADVECTION 

In we showed that advection of a large-scale mag- 
netic field will dominate over diffusion in the nonturbu- 
lent surface layer of an accretion disk and that, if equa- 
tion ([U is satisfied, the advection can lead to a concen- 
tration of magnetic flux in the inner region of the disk. 

In this section, we present a more rigorous analysis of 
the conditions under which equation ^ is satisfied. In 
we consider the forces that act between the main 
body of the disk and the material at the base of the 
nonturbulent region, which determine whether or not this 
region is actually a true "surface layer" that participates 
in some way in the accretion fiow (i.e., whether or not 
Vrh < there; clearly, at some height above the disk, 
the radial velocity may no longer be inward). In !j4.21 
we use these results to derive conditions on the magnetic 
field geometry that allow inward advection of magnetic 
fields to proceed, and in M.31 we derive more general 
constraints on the field strength and ionization fraction 
that are required for a highly conducting nonturbulent 
surface layer to exist in the first place. (In general, the 
above constraints are weak, and whether or not advection 
can occur is primarily a question of geometry; a vertical 
or "dipole-type" seed field provided in the outer regions 
of the disk will often be sufficient to advect inward on its 
own.) Finally, in t j4.41 we consider the ultimate outcome 
of the advection of a large-scale magnetic field and the 
types of disks in which it might generally be expected to 
occur. Readers not interested in the detailed analysis we 
present in the rest of this section may wish to skip to the 
conclusions of this paper in fj5l 
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4.1. Forces Acting on the Nonturbulent Surface Layer 
of the Disk 

For an axisymmetric accretion disk in which the 
specific angular momentum profile is time-independent 
(which can be true even in a time-dependent disk if or- 
bits are nearly circular), a general equation for the local 
radial velocity is 



_dr 



r4>) 



rp- 



d_ 

dz 
"5 



(r T^z) + r pvz-^ 



(5) 



This expression is obtained by combining the conser- 
vation of mass and conservation of angular momentum 
equations for a magnetized fluid. Here, p is the mass 
density and T^^ and T^z are components of the stress 
tensor, including both large-scale magnetic and small- 
scale turbulent stresses. If we evaluate this equation at 
z = h (the height in the disk where turbulence is first 
suppressed) and if the gravitational and centrifugal forces 
are assumed to balance in this region (i.e., if we assume 
circular orbits in Newtonian gravity), then to first order 
in h/r, 



Vrh 



d_ 

dr 



[r'BrBA, 



d_ 

dz 



2TTrphVK 



1 Vzh, (6) 



where is the Keplerian velocity on the equatorial 
plane. The first term represents the effect of the stress 
due to the time-averaged magnetic field, and the second 
term represents a centrifugal effect that drives outflowing 
material away from the inner disk. Competition between 
these two processes determines the vertical profile of the 
radial velocity in the nonturbulent region of the disk. 

The analysis leading up to equation ^ is quite general, 
and we therefore make liberal use of it in the following 
sections. The only exception to its generality is the as- 
sumption of circular orbits, which may not be valid if the 
base of the nonturbulent region occurs at a height h in the 
disk where the density is so low that radiation pressure 
or magnetic forces begin to become important in the ra- 
dial momentum equation. This is effectively a constraint 
on the magnetic field strength, and we therefore discuss 
it further in ! j4.3l In the meantime, note that equation 
([5|) does not assume circular orbits (provided that the 
disk is stationary) and that the derivation leading from 
equation ([5]) to ([B]) will still be roughly valid provided 
only that the spatial derivatives of the azimuthal veloc- 
ity Vfj,h in the nonturbulent region are of the same order 
of magnitude as their Keplerian counterparts. 

4.2. Conditions on the Magnetic Field Geometry 

A straightforward way to derive a condition on the 
magnetic field geometry is to combine equations ([4]) and 
([6]). If we do this and assume that the stress due to 
the time-averaged magnetic field removes angular mo- 
mentum from the disk surface (i.e., attempts to drag the 
surface layer inward with the main body of the disk) 
in any amount, then a sufficient condition for growth of 
magnetic flux in the inner disk is 



Brh 

Bzh 



-] VzH > 0, 

r 



(7) 



provided that magnetic field with the appropriate polar- 
ity is available to accrete. 

When equation ([7|) is satisfied, magnetic flux growth 
can occur through a combination of radial and vertical 
advection at the surface of the disk. However, we are 
primarily interested in radial advection, which is the only 
sustainable way in which magnetic field lines anchored in 
the surface layer can build up flux in the inner disk. In 
particular, when equation ^ is satisfied, we can identify 
three regimes of interest: 

1. If Vrh < < [Brh/ Bzh) Vzh, magnetic flux growth 
occurs through a combination of inward radial ad- 
vection and vertical advection at the surface of the 
disk. 

2. If Vrh < {Brh/ Bzh) Vzh < 0, magnetic flux growth 
occurs via inward radial advection, even though it 
is partially opposed by vertical advection. 

3. If < Vrh < {Brh/ Bzh) Vzh, magnetic field is ad- 
vected radially outward at the surface of the disk, 
but magnetic flux growth still occurs in the inner 
disk because of vertical advection at z « ft,. 

We are primarily interested in the first two regimes, 
where magnetic flux growth occurs at least partially due 
to inward radial advection, and where equation ([7]) does 
not necessarily apply. Thus, for the rest of this section, 
we ignore the third regime and derive more stringent con- 
ditions that specifically guarantee inward radial advec- 
tion. 

As discussed in fj2l the first regime is likely to be 
more relevant than the second (jUstvugova et al.l 119991 
[2000I) . and we therefore consider it now, returning to 
the second regime at the end of this section. We thus 
have {Brh/ Bzh) Vzh > and require Vrh < for in- 
ward radial advection. If we use equation ^ and define 

Hb = [d\n{B^hBzh) /dz]~ as the scale height of the 
vertical magnetic stress, we find that a sufficient condi- 
tion is 



B^hBzh ^ - 



3 Vzh Hb 

2 (Vr) H 



Ph_h 
Po r 



K 



(8) 



where po is the mass density on the equatorial plane, 
M is the local mass accretion rate, ri/j- = v^^/r is the 
Keplerian angular velocity, and {vr) = M / {AnrpoH) is 
an appropriately height-averaged inward radial velocity 
in the main body of the disk (i.e., the standard radial 
velocity in a one-dimensional vertically-integrated disk 
model). In interpreting this equation, it is instructive to 

note that {MflK/r)^^^ « 1-8 x 10^ m-^^^m^/'^f-^^'^ G is 
a fiducial field strength, but a maximum one, since h/r 
and especially ph/po can be very small parameters.^ 

Equation ([5]) is a "sufficient" condition for magnetic 
field advection in the sense that it makes the conserva- 
tive assumption that large-scale magnetic stresses in the 

^ Although we can generally assume h ^ H within a factor of 
a few, the distinction between those two heights must be retained 
when evaluating the mass density, and it is important to use the 
correct value pj, which appears in equation ((HJl. This is because 
the mass density typically falls off very sharply with height, and 
thus may be many orders of magnitude smaller than both pu 
and Po- 
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vertical direction are the only way in which the surface 
layer can be dragged inward. In particular, it does not 
include the large-scale BrB^ stress at the disk surface, 
which transports angular momentum radially and also 
tends to give Vrh < (although for a thin disk, the ver- 
tical stress is usually most important). 

The apparent dependence on the field strength in equa- 
tion ([8]) disappears when the equation is analyzed more 
carefully. In particular, we show in Appendix |A] that 
the definition of h as the height in the disk where MRI 
turbulence is first suppressed allows us to write 

^ - (9) 

where Bzh ~ \Bzh\ + [H/r) \B^h\ represents a magnetic 
field strength that is roughly equal to the time-averaged 
vertical field Bzh^ and we define po = Po {HVLk)^] for a 
weak field, po is roughly equal to the thermal pressure 
on the equatorial plane of the disk. Combining equations 
(HI) and ® and assuming /i ~ iJ, we find that a sufficient 
condition for inward advection of magnetic fields is 



BthhB 



> 



Hb Vzh 



(10) 



which has no direct dependence on the field strength; 
as long as the geometry is favorable, arbitrarily weak 
magnetic fields can provide enough stress to drive inward 
radial advection at the surface layer of the disk, and the 
fields will therefore be advected inward and compressed. 
The essential physical point is simply that the magnetic 
field must be strong compared to the gas at z ~ /i (in 
order to suppress the MRI), so it is therefore able to drive 
accretion at this location, regardless of how weak it is in 
an absolute sense.'' 

If we ignore the sign of B^hBzh, equation (fTO)l is rel- 
atively trivial to satisfy. For example, we can estimate 
the ratio of vertical magnetic stress to energy density 
on the left hand side of equation ([TI?)) that might arise 
naturally in a disk (i.e., without an externally-imposed 
seed field) by looking at numerica l simulations of th e 
MRI. We focus on the work of ,Miller fc Stond (|200ClD . 
who studied a vertically-stratified disk in the shearing 
box approximation and tabulated the properties of the 
magnetic field in the nonturbulent corona above the 
disk. We find typical values of \B^hBzh\ /Bl,^ > 0.05 
in these simulations. This should be compared to the 
right hand side of equation (fTO|) . whose magnitude is 

given by ~ a {H/r)^ [Hb/H) \vzh \ / (vr) when the large- 
scale field is dynamically weak (or, alternatively, Hb/t 
times the ratio of Vzh to the disk sound speed). This 
is clearly a very small number if we make the approxi- 
mation that Hb ^ H (i.e., that the scale height of the 

This statement is independent even of the Spitzer difEusivity; 
our expression for v^h in equation ^ means that we can rewrite 
the condition Rem 3> 1 (which must be satisfied at the base of 
the nonturbulent region in order for advection to dominate the 
Spitzer diffusivity there) as \(H/Hb) B^hB;^h/Bl^ + v^^/^kI > 
10-"m-5/8m-3/8f5/8/7^/* (l02_H"/r)"^. Clearly, this equation 
will be satisfied in almost any accretion disk provided that equation 
moi l is not pathologically close to an equality. More important 
for our purposes, there is no dependence on the magnetic field 
strength; even when microscopic effects such as Spitzer difEusivity 
are taken into account, arbitrarily weak fields appear capable of 
advecting inw ard a long the surface layer of a fully ionized accretion 
disk (but see 



twisted toroidal field can be comparable to or smaller 
than that of the mass density); the validity of this ap- 
proximation is discussed in Appendix |B] Intuitively, the 
approximation Hb ^ H may be thought of as arising 
from the presence of a voltage source (Keplerian shear) 
that is applied in the radial direction, with the result- 
ing current confined to flow in a region above ~ h (the 
height at which the plasma first becomes highly conduc- 
tive) but below ~ a iewxH (the height at which the mass 
density becomes low enough so that orbits are no longer 
circular and therefore the applied voltage is significantly 
reduced). 

The ease with which equation pO|) can be satisfied sug- 
gests that not only will we have Vrh < in a typical 
accretion disk, but the accretion fiow may also reach a 
steady state in which the main, turbulent body of the 
disk drags the nonturbulent surface layer inward at the 
same speed as itself; i.e., \vrh\ ~ {vr}- In fact, if we start 
with equation ([6|) and go through the same analysis as 
above but require Vrh ^ — (vr) rather than Vrh < 0, then 
instead of equation pI7|) we obtain 



B^hBzh ^ Hb {vr) 



zh 



H Hil 



K 



1 



3/i Vz 



(11) 



where the right hand side is typically ~ uHb/t for a 
dynamically weak field. Like equation (|10|). this condi- 
tion is modest assuming the disk is thin, and thus we 
may expect that advection of magnetic fields proceeds at 
the same speed as turbulent accretion in the main body 
of the disk. In fact, equation pT|) suggests that advec- 
tion of magnetic fields in the surface layer could proceed 
faster than the disk accretion speed, but as we discuss in 
Appendix [b1 this is unlikely to be sustainable. 

The only qualification to what we have said so far con- 
cerns the sign of B^^Bzh- In particular, equation ([8|) 
shows that B^^Bzh ^ is required in order for the field 
to advect inward; the exact condition is d {B^Bz)^^ /dz < 
0, which states that the large-scale magnetic field must 
remove angular momentum from the nonturbulent sur- 
face layer. This is a strict requirement for disks in which 
Vzh > 0. We expect that d {B^Bz)i^ /dz < will be satis- 
fied in many parts of an accretion disk, but not necessar- 
ily all. It will be satisfied in regions where the magnetic 
field has a dipole-type symmetry (with Br and B^ odd 
functions of z and Bz an even function) , which is often as- 
sumed for the large-scale magnetic field advected inward 
in an accretion disk (see Figure [1]). However, in a region 
of the disk with a quadrupole-type field symmetry (with 
Br and B^ even functions of z and Bz an odd function) , 
as may occur when the large-scale field extending out 
of the disk is ge nerated primarily by mag netorotational 
turbulence (e.g.. iBrandenburg et al.lll995f ), some regions 
will hkely have d{B^Bz)f^ /dz > 0. In these regions, an- 
gular momentum will not be removed vertically from the 
surface of the disk, and inward radial advection of the 
magnetic field may be difficult to sustain. Correspond- 
ingly, w e find that one of th e simulations described in 
detail in iMiller fc Stond (|2000) has d{B^Bz)f^/dz < 0, 
but the other does not, so that the material (and mag- 
netic field) in the nonturbulent region may not advect 
inward in this second case. 

On the other hand, if there is a weak vertical seed 
field in the outer part of the disk, shear and MRI tur- 
bulence will create azimuthal field from it, and the con- 
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dition d {B^Bz)^^ / dz < should be satisfied automat- 
ically, while the condition in equation (jf Op will be un- 
changed. We therefore view this as the most favorable 
way to induce magnetic field advection in a disk. The 
large-scale magnetic field introduced into the disk at 
large distance may be supposed to come from the inter- 
stellar medium (in the case of supermassive black holes) 
or from a companion star (in the case of X-ray binaries 
and other stellar-mass systems). 

To conclude this section, we consider the second regime 
for magnetic field advection alluded to in our earlier dis- 
cussion of equation in which the geometry is such 
that {Brh/ Bzh) Vzh < 0, and inward radial advection re- 
quires Vrh < {Brhl Bzh) ^zh in ordcr to overcome vertical 
advection and produce a concentration of magnetic flux 
in the inner region of the disk. If the disk is thin, with 
3/i/r <C \Brh/ Bzh\, a sufficient condition for this to oc- 
cur can be obtained by replacing h/r with —Brh/SBzh in 
equation ([5]) and propagating this change through subse- 
quent expressions. In particular, equation (|10p becomes 

BijjhBzh y Hb f Brh\ Vzh 

where we have {—Brh/ Bzh) Vzh > by definition so that 
the right hand side of the equation is, in general, posi- 
tive. This condition can be satisfied for \Brh/Bzh\ ^ 1 
provided that Vzh is not too close to the sound speed. 

4.3. Conditions on the Magnetic Field Strength and 
Ionization Fraction 

From the analysis of the previous subsection, an arbi- 
trarily weak seed field threading an accretion disk should 
be able to advect inward along the disk surface. But what 
really happens for arbitrarily weak fields? Are there field 
strengths below which some of our underlying assump- 
tions in this paper break down? 

An important assumption in this paper is that the re- 
gion above the disk is nonturbulent and, therefore, highly 
conducting; it is this region in which the magnetic field 
can advect inward. Clearly, a nonturbulent region is 
likely to exist somewhere above an accretion disk, but 
the question is whether it occurs at a low enough height 
to be treated as the "surface layer" of the disk, as we do 
in this paper. Equation ([9]) and the usual assumption 
that the mass density decreases with height much more 
rapidly than the magnetic energy density suggests that 
the nonturbulent region should occur within a few scale 
heights, even for a very weak seed field on the equatorial 
plane. However, if the magnetic energy density begins 
to drop off rapidly with height, the turbulence may not 
be suppressed until a very large distance above the disk. 
This may be what happens in the radiation-dominated 
simulations of lTurnerl (|2004D . where the magnetic energy 
density begins to fall off at z > 3H, and there is no clear 
evidence for a nonturbulent region anywhere within the 
simulation domain (which extends out to z ~ 8H). 

It is difficult to predict when this type of behavior will 
occur, but when it does, our assumption that the nontur- 
bulent region occurs "within the disk" may break down. 
In particular, orbits may not be circular, so that equa- 
tion ([6]) is no longer strictly valid. Considering radial 
force balance and using equation ([5]), we find that if the 
disk is sufficiently thin, magnetic forces are unlikely to 



be strong enough to disrupt circular orbits at the base of 
the nonturbulent region (although they may certainly do 
so higher up in the corona); the condition for magnetic 

forces to be negligible at z « /i is Bzh/B^ {H/r)^ 
and Blf^l \Brh Bzh\ ^ H/r, both o f which are easily 
satisfied by the iMiller fc Stond ()2000l ) simulations (here 
Blf^ ~ Bly^ was defined in ^A.2\ and in the second expres- 
sion we have assumed that the scale height of \BrhBzh\ 
can be approximated as H). Radiation pressure is 
therefore the only realistic concern; in order for orbits 
to remain circular, we require Ph^j^ to dominate over 
the radiation pressure. Assuming the temperature at 
z ~ h \s given by the effec tive surface temperature of a 
IShakura fc SunvaevI (|1973[ ) disk, we can use equation ([9]) 
to derive 

Bzh » 10-2 m-^'^m^'^ff^fl'^ (iQ^H/r) G (13) 

as the condition for circular orbits, where = f/10'^ and 
mg = m/10* (i.e, we have scaled the fiducial value to that 
which would occur at a distance of 10'^ Schwarzschild 
radii from a supermassive black hole of mass IO^Mq). 

Limits on the magnetic field strength resulting from 
equation (fT3l) are plotte d in Figure [3] for typical 
IShakura fc SunvaevI (|1973f l accretion disks with a « 
10"*, which we take to be a worst-case lower limit for the 
turbulent stress (note in any case that the dependence on 
a is very weak). In fact, a is not independent of the field 
strength; using the relation a ^ O.SSq/Sttpo between 
turbulent magnetic stress and turbulent magnetic energy 
density Bq /Stt in the main body of th e disk (which is a ro- 
bust result of MR I simulatio ns: e.g.JHawlev et a,l lll995t 
ISano et al.l |2004 [Blackma nTPenna. fc Varnierd \20m . 
we can rewrite equation (|13p as a constraint on the field 
geometry: 

^»AxlO-^(^P^) 'm'^mf-^U (14) 
Bq Vgcm-V 

Limits on this ratio for typical IShakura fc SunvaevI 
(|1973f ) accretion disks are plotted in Figure [H again as- 
suming a « 10^'*. In the radiation-dominated region of 
the disk, po oc a~^, so the right hand side of equation 
(fH|) has no dependence on the magnetic field strength, 
whereas in the gas pressure-dominated region, the de- 
pendence occurs through a~^^^^ , which we have fixed to 
the assumed worst-case value. Thus, although equation 
(jl4p technically represents a joint limit on the equatorial 
plane field strength and the vertical magnetic geometry, 
the condition on the geometry is more important. 

In a conservative analysis, equation (|13p can be com- 
pared to the large-scale magnetic field that might be sup- 
plied to the outer disk (from, say, the interstellar medium 
or a companion star) in order to determine whether or 
not orbits are circular in this region. However, equation 
(jl4p may be a more appropriate expression if we allow 
for the possibility that buoyant rising of magnetic fields 
from the turbulent disk into the nonturbulent corona can 
a ffect the magnitude oi Bzh for example, the simulations 
of lMiller fc Stmia (|200nh have BjjB^ > 0.02 (even with- 
out any net imposed vertical field) and thus should easily 
satisfy equation ([H]) . 

As mentioned in ! j4.11 even when orbits are not circu- 
lar, equation ^ and the subsequent analysis may still 
be valid to an order of magnitude, and so we do not view 



8 



Rothstein & Lovelace 



II 10^ 

■dS 10' 
■-2 

o;^ 10° 

i| 10-^ 
2 10-^ 













- 






- 





















10^ 



10^ 



10' 



10" 



10' 



lO*" 



Radius (Schwarzschild radii) 



Fig. 3. — Minimum values of the vertical magnetic field B^^ that 
are required in order for the large-scale field to overcome radiation 
pressure and be advected inward in the nonturbulent surface layer 
of an accretion disk. Constraints are shown for black hole masses 
of lOMp) and lO^Mp, and dimensionless accretion rates between 



rh = 0.1 and m = 10 (note that the Eddington accretion rate 
corresponds to rh = e~^, where e is the radiative efficiency of 
the disk and e fa 0.06 for a simplified treatm ent of accretion onto 
a Schwarzschild black hole). We assume a IShakura &: SunvaevI 
l|1973i ) disk with /« = 1 at the inner boundary and a = 10~* 
as a worst-case lower li mit for the turbulent stress. Solid lines 
show regions in which the IShakura fc Sunvaev I l[l973) assumptions 
about pressure and opacity sources in the disk are good to at least 
50%; dotted lines are used to connect through transition regions 
in which the IShakura fc SunvaevI I I1973I) solution breaks down. 



the constraints on the magnetic field discussed here as 
fundamental lower limits. On the other hand, one might 
reasonably expect that the effect of radiation pressure is 
to drive material radially outward from the hot inner disk 
and thereby prevent the nonturbulent surface layer and 
its associated magnetic field from accreting; in that case, 
circular orbits can indeed be viewed as a strict require- 
ment. Conversely, other mechanisms that might produce 
noncircular orbits, such as advectio n-dominated flows 
at cither high or low accretion rates ijAbramowicz et all 
[1988: Na ravan fc Yi 1994, 199S), tend to enhance the in- 
ward radial velocity and therefore do not prevent inward 
advection of magnetic fields from occurring. In these 
cases, equation ([TU)) should have vk replaced by the ac- 
tual azimuthal velocity v^^ which makes the equation 
more difhcult to satisfy but does not change our overall 
conclusions. 

Even if orbits in the nonturbulent region are circular, 
other constraints may come about that could affect the 
ability of this layer to advect magnetic fields inward. In 
particular, our discussion so far has implicitly assumed 
that the disk is fully ionized. We do not discuss partially 
ionized disks in depth in this paper, but we point out two 
important issues. First, the diffusivity associated with 
electrons scattering off of neutrals (as well as other non- 
ideal MHD effects) may become more important than the 
electron-ion Spitzer diffusivity, and second, there may 
not be enough free electrons in the nonturbulent surface 
layer of the disk to support the advected magnetic field. 

We cons ider non- ideal MHD effects first. Following 
iPandev fc' Wardlc ( 2007) , a more general version of equa- 
tion (HI) which includes the relevant non-ideal MHD 
terms is 



dB 

'dt 



= V X 



v'xB-i^J 

c 



pI (JxB)xB 



P PiVinC 



(15) 
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Fig. 4. — Same as Figure [3| but for the minimum values of the 
ratio B'^i^/Bq (between the energy density of vertical magnetic 
fields at the disk surface and turbulent fields on the equatorial 
plane) that are required in order for advection of the large-scale 
magnetic field to overcome radiation pressure. 

where v' = v — 3 /upC is the plasma velocity modified by 
the Hall drift, J = [c/Att) V x B is the current density, 
Ue is the electron number density, e is the proton charge, 
Pn and Pi are the neutral and ion mass densities, and Vm 
is the ion-neutral collision frequency. The last two terms 
in equation represent the effects of Ohmic diffusion 
(i.e., the scattering of electrons off of ions and neutrals) 

and ambipolar diffusion, r espectively. 

We use equations in iDraine. Roberge. &: Dalgarnol 
(|1983f l and iBalbus fc TerquemI (j2001[ ) to estimate the 
importance of the non-ideal MHD effects. For simplic- 
ity, we assume a priori that the plasma is weakly ion- 
ized and that the electron and ion number densities 
are the same. We take "worst-case" approximations 
of quantities involving the magnetic field (for example, 
|V X B|^ ^ Bh/H), assume that the plasma tempera- 
ture is gi ven by the effect ive surface temperature of a 
[Shakura fc SunvaevI (|1973f l disk, and make use of equa- 
tion © to substitute for the density when needed. Fidu- 
cial values are given for a distance of 10'^ Schwarzschild 
radii from a IO^Mq black hole, where the surface tem- 
perature is of order 2,000 K and thus the disk should 
be weakly ionized. For diffusive processes, we can calcu- 
late an effective magnetic Reynolds number Rem in the 
nonturbulent surface layer of the disk (see ^J2|). When 
Re„i 3> 1, advection of the magnetic field will dominate. 
For Ohmic diffusion, we find that this condition can be 
approximately written as 



> 5 X 10^ 



-15 



-9/8 — 7/8 



Ohm 1 



(16) 



where rii/n„ is the ionization fraction at the base of 
the nonturbulent region (i.e., the ratio of ion to neu- 
tral number densities) and we have defined a parameter 

Dohrn = (10a)"^mi/8/*^^ (lO^if/r)"^^^,-! which con- 
tains terms of lesser importance. Here, Us = \vrh \ / (vr) 
is likely to be ~ 1 as per our discussion in ii4.2l For 
ambipolar diffusion, the condition becomes 

2 

> 10"^ 



G 



-lj.-5/2^ 

3 ^ambi 



(17) 



where _D„ 



= (10«)'' {10^ H/r) p' \Bh/5Bj'Ur' 
and p' = /i (1 + nin/mi); here, p is the mean mass per 
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particle expressed in units of the proton mass and m„/mi 
is the ratio of neutral to ion masses. 

We therefore see that ambipolar diffusion imposes a 
joint constraint on the ionization fraction and magnetic 
field strength that is generally much more important 
than the limit imposed by Ohmic diffusion; this is a di- 
rect result of the low densities expected in the nontur- 
bulent surface layer that we consider in this paper. For 
a partially ionized disk that is able to overcome Ohmic 
diffusion and satisfy equation (jl6p . magnetic fields can 
advect inward, but equation (jl7p shows that stronger 
initial seed fields than those indicated in equation (|13p 
and Figure [3] (for the constraint imposed by circular or- 
bits) are generally required. However, in these applica- 
tions it is important to keep in mind that the ionization 
fraction that appears in the above equations is calcu- 
lated at the surface of the disk (as opposed to the main 
disk body), where interstellar cosmic rays may inflate the 
ionizati on rate above the values that one would usually 
expect ()Gammielll996[) . 

From equation (fT5|) . we see that the Hall effect, which 
occurs due to the drift of the magnetic field with respect 
to the ions as it is carried along by the electrons, may also 
be important, not only for partially ionized disks but also 
for fully ionized ones. However, the Hall effect does not 
cause diffusion of the magnetic field and therefore does 
not necessarily oppose inward advection. In fact, if the 
radial component of the current density is positive at the 
surface of the disk. Hall drift enhances inward advection 
rather than opposes it; this occurs when dB^h/dz < 0, a 
condition which is already met by the dipole-typc fields 
discussed in ^14.21 In geometries where the Hall effect 
opposes inward advection, however, the condition for it 
to be negligible is 



rin. 



B 



^ ) > 2 X 10-^ m« 



-3/2 



Dnaii, 



(18) 



where Dnaii = {Wa)'' U (lO^iJ/r) \^l \Bh/bB,h\ U-\ 
This equation also applies for a fully ionized disk if we 
replace ni/un by « 1/2. Note that the Hall effect is 
generally more important in a fully ionized accretion disk 
than the electron-ion Spitzer diffusivity discussed in Sj2l 
and in some cases (especially for stellar mass objects), it 
may even impose a stronger constraint on the magnetic 
field strength than that imposed by equation (|13p for 
circular orbits. 

We finally consider the possibility that in the diffuse 
surface layer of a weakly ionized disk, there may not be 
enough free electrons available to carry the current that 
is necessary for supporting the large-scale magnetic field. 
Making similar assumptions as above and also assuming 
that the velocity of electrons is limited to their thermal 
speed, we find that 



rin 



G 



> 9 X 10 



-12 



zh\ 



7/8 . 1/8.13/8 .1/8 



(19) 



is required to produce enough current to support the 
magnetic field, which is generally a weaker constraint 
than that imposed by ambipolar diffusion. (Note that 
stronger magnetic fields are actually easier for the non- 
turbulent surface layer to support; this is because strong 
fields suppress turbulence at a lower height above the 
disk, where the electron density is larger.) 



4.4. The Outcome of Advection 

In summary, we find that weak, large-scale magnetic 
fields can be advected inward in the surface layer of an 
accretion disk. The most important condition necessary 
for this is d {B^Bz)i^ jdz < (i.e., the vertical magnetic 
stress must extract angular momentum from the disk sur- 
face). This may occur most easily if a weak, dipole-type 
seed field is supplied at the outer regions of the disk. 
This field threads the disk and acts as a catalyst for pro- 
duction of the toroidal fields that in turn provide the 
required geometry for further inward advection. 

Although the constraints on the field strength dis- 
cussed in ij4.3l are modest, we can nonetheless imagine 
that a weak, dipole-type magnetic field (relatively close 
to the lower limit) is supplied in the outer region of the 
disk and ask the following question: How easy will it be 
for the field to remain at the required level as it is ad- 
vected inward and compressed? From equation (jl3p . we 
find that roughly, the vertical field must increase faster 
tha n Bz oc r~'^/^ in the gas pr essure-dominated region 
of a IShakura fc SunvaevI (|1973D disk in order to satisfy 
the requirement for advection to continue; this limit be- 
comes Bz cx r~^/^ in the inner, radiation-dominated re- 
gion (where H is approximately constant). This growth 
seems difficult to sustain, which suggests the interest- 
ing possibility that an advected magnetic field may tem- 
porarily "stall" at some large radius when it becomes too 
weak to advect inward, and advection will only continue 
once enough magnetic fiux has built up at this location. 
On the other hand, as we have seen in discussing equa- 
tion (fTi|) . even if the advected field does not increase in 
strength fast enough on its own, it may be reasonable 
to assume that the local dynamo can produce strong 
enough surface fields (through magnetic buoyancy) to 
meet these conditions, and thus advection will continue 
in either case. 

As discussed in 21 a sufficiently strong magnetic field 
threading the disk can lead to inward radial advection 
of the field driven by the main body of the disk too, 
due to extraction of angular mom entum from the mai n 
body of the disk to a wind or jet (jLovelace et al.|[l994f ). 
However, the field strengths required for this process are 
much larger; the field must drive accretion in the dense 
body of the disk (rather than in a low-density surface 
layer). In fact, advection of magnetic fields in the surface 
layer is always more efficient, provided that the disk is 
thin and the vertical field is subequipartition. 

Although advection in the main body of the disk must 
take place fast enough to overcome turbulent diffusion, 
advection in the nonturbulent region is limited by the 
much smaller speeds associated with the microscopic dif- 
fusivity. Magnetic fields can advect inward in the non- 
turbulent surface layer even when \vrfi\ <C (vr), and the 
advection can therefore take place on timescales longer 
than a viscous timescale. This will lead to a gradual 
buildup of field in the inner disk. In the case of an- 
gular momentum extraction from the main body of the 
disk, advection of magnetic fields always occurs at least 
as fast as the viscous timescale, and sustained advec- 
tion requires that the magnetic field be strong enough to 
completely overwhelm turbulent diffusion; otherwise, a 
steady state will be reached in which only a small con- 
centration of magnetic field develops in the inner disk. 
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Thus, we may expect a disk to experience gradual ad- 
vection of large-scale magnetic fields on long timescales 
in the surface layer, wi th occasional brief bu rsts of fast 
or implosive accretion (|Lovelace et al.|[l994l) associated 
with the presence of a strong magnetic field in a local 
region of the disk that extracts angular momentum from 
the main disk body. 

It is interesting to note that advection of magnetic 
fields in the surface layer is more efficient for thin disks 
than thick ones, as can be seen, for example, from equa- 
tion ([8]). This is because the surface layer is tightly cou- 
pled to the main body of the disk rather than being part 
of a detached corona. On the other hand, turbulent ad- 
vection of magnetic fields in the main body of the disk 
is more efficient if the disk is thick (since Rem ~ H/r, 
as discussed in Sjjl- We therefore conclude that mag- 
netic field advection is possible for many disks but most 
efficient for thin ones. 

5. CONCLUSIONS 

This paper reanalyzes the advection of a large-scale, 
weak magnetic field in an accretion disk. We consider the 
vertical structure of the disk, which strongly influences 
the vertical profile of the conductivi ty, as pointed out by 
iBisnovatvi-Kogan fc Lovelacd ()2007[ ) . In the thin, diffuse 
surface layers of the disk, the magnetic energy density is 
large enough compared to the thermal energy density 
that magnetorotational turbulence is suppressed. As a 
consequence, magnetic field lines threading the surface 
layer can be advected inward with the main body of the 
disk, without being opposed by turbulent diffusion. 

No special conditions are required for the field to be ad- 
vected inward except that it meet the rather modest con- 
straints in 21 The required field strengths are relatively 
weak, and the primary constraint on the field geometry 
is simply that it must help the nonturbulent surface layer 
accrete inward (i.e., the vertical magnetic stress must ex- 
tract angular momentum from this layer, in virtually any 
amount). This can be accomplished either via coupling 
between the main, turbulent body of the disk and the 
surface, or via a wind or jet. The simplest way for this 
condition to be met is if the accretion flow is provided 
with a weak, large-scale vertical seed field threading the 
outer region of the disk (which could come from the inter- 
stellar medium or a companion star), although in some 
cases, the proper geometry may be attained entirely as a 
result of fields produced via the local magnetorotational 
dynamo. Once a weak, large-scale field with the proper 
geometry is in place, the field will be advected inward 
along the disk's surface layer and strengthened as it is 
compressed along with the accretion fiow. 

The presence of a large-scale magnetic field anchored 
in the surface layer will drive strong magnetorotational 
turbulence in the main body of the disk, which can con- 
sequently produce values of the a parameter (i.e., the 
turbulent stress) large enough to match observational 
constraints. We find that typical vertical fields on the 
order of a few percent of equipartition are required for 
this to occur. Because the field is advected inward and 
anchored in the surface layer, there is no need to worry 
about maintaining the required vertical fields via in- 
ternal magnetorot ational fluctuations (as suggested by 
iPessah et al.l[200l . We propose that long-term changes 
in a (in response to the history of magnetic fleld ad- 



vection) should be explored as a possible source of the 
long-term evolution in the variability patterns seen in 
the light curves of X-ray b inaries such as GRS 19154-105 
(see, e.g.. lTagger et al.| ]2004V 

The mechanisms we discuss in this paper are relevant 
to many different kinds of accreting objects. In the outer 
part of the disk (far away from the central star or black 
hole), our work should be applicable provided only that 
the disk is sufficiently ionized (see H4.3p . Closer to the 
inner part of the disk, our work is most obviously ap- 
plicable to black holes, where the large-scale magnetic 
field arises entirely within the accreting plasma. How- 
ever, it may also be relevant to the case where a large- 
scale field from a magnetized central star penetrates 
the disk; previous suggestions that outward radial dif- 
fusion of the field may b e an important process i n such 
systems fe.g.. [ Lovelace. R omanova. & Bisnovatvi- Kogan 
1995*:'Bar dou fc Hevvaer ts 1996; Aga pitou fc Papaloizou 
2000; Uzde nskv. Konigl. fc Litwin 20(^ should be exam- 
ined in light of our work. Finally, the applicability of our 
work to radiation-dominated regions of the disk may re- 
quire further analysis, because it is unclear whether the 
vertical structure of these regions can support the inward 
advection of magnetic fields (although we argue here that 
it can). Future numerical simulations may eventually be 
able to address this point. 

Also, we have assumed axisymmetric disks in this pa- 
per, but the equations we have derived apply more gen- 
erally if they are averaged over azimuth (with the in- 
troduction of appropriate correction factors). Thus, a 
magnetic field with an azimuthally- averaged strength and 
azimuthally- averaged geometry that roughly meets the 
criteria in f}4]is likely to be able to advect inward in an ac- 
cretion disk. Nonaxisymmetric adve ction of a large-scale 
magnetic field (e.g., as envisioned bv lSpruit fc Uzdenskv! 
[2005) is clearly consistent with the mechanism we have 
proposed in this paper, but we have also shown that 
advection is equally plausible in an axisymmetric disk, 
without having to assume any special conditions. More 
numerical simulations that investigate disks in which a 
large-scale magnetic field is supplied at the outer bound- 
ary are clearly needed. 

Finally, our work has important implications for mod- 
els of jet formation which require strong, large-scale mag- 
netic fields to exist over a region of the inner accretion 
disk. As suggested more than 30 years ago, these mag- 
netic fields can arise in the inner disk via a very simple 
process: advection of a weak field from outside. This 
opens up the possibility for magnetically-dominated out- 
flows (i.e., Poynting jets) to exist in the inner regions 
of disks around a wide variety of accreting objects. In 
addition, the radial stretching of fleld lines produced by 
advection may al low winds accelerated by the magneto- 
centrifugal effect (jBlandford fc Pavii3ll982f ) to exist over 
a wide range of radii. This is in contrast to the results 
of numerical simulations in which no magnetic fields are 
supplied at the outer boundary, and jets only form in 
an extremely narrow inner region where the energetics 
may be dominated by relativistic effects that require the 
presence of a rotating black hole. 

ACKNOWLEDGEMENTS 

We thank G. S. Bisnovatyi-Kogan, M. M. Ro- 
manova, and I. G. Igumenshchev for valuable discus- 



Advection of Magnetic Fields in Accretion Disks 



11 



sions. D. M. R. is supported by an NSF Astronomy by NASA grants NAG5-13220 and NAG5-13060 and by 
and Astrophysics Postdoctoral Fellowship under award NSF grant AST-0507760. 
AST-0602259. The work of R. L. was supported in part 

APPENDIX 

A. PHYSICAL CONDITIONS AT THE BASE OF THE NONTURBULENT REGION ABOVE AN ACCRETION DISK 

In this section we derive equation © , which defines the base of the nonturbulent surface layer of the accretion disk 
and is used throughout the main body of the paper. 

We are interested in conditions at z = /i, the height in the disk where the magnetorotational instability (MRI) 
is first suppressed. The fundamental condition for suppression of the MRI is that the Alfven speed must be large 
enough so that fluid elements linked by the magnetic field will be drawn back together faster than orbital shear drives 
them apart; in a WKB (small wavelength) analysis, M RI modes with wavenumb er k are found to be suppressed when 
k • va > ^K, where va is the local Alfven speed (e.g., iBalbus fc IIawlevlll998f) and we have assumed circular orbits 
(relaxing this assumption would simply replace fix with the actual local angular velocity of the disk). If we take the 
WKB approximation to its limit and make t he usual assumption tha t the MRI will be completely suppressed when no 
unstable wavelength fits within the disk (e.g.. lBalbus &: IIawleylll99"ll ) and further consider that the largest wavelengths 
that fit within the disk are of order (a few) x H in the vertical direction and ~ (a few) x r in the azimuthal direction, 
then applying this condition to the material at z = /i shows that MRI turbulence will be suppressed when 

{vA.)h + - ivA4,)h > HQk, (Al) 



where (waz)/; = B / \/ Ph and {va^,)^^ is defined equivalently. Rearranging this equation gives, approximately. 



where the equality holds at z « /i (the height at which turbulence is first suppressed), but the equation also applies 
more generally if h is redefined to be any height within the nonturbulent region. Here, we have defined a modified 
vertical magnetic field strength Bzh = \Bzh\ + {H/i') \B^h\- In the main body of the paper, we generally assume 
Bzh ~ \Bzh\i but the full expression should be used when the toroidal field is extremely strong {Bzh ~ \Bzh\ appears 
to be a good approximation for the magnetic fields seen in the simulations of Miller fc Stone. 200Q . however). 

If we define po = Po (H^k) , equation (jA2p becomes 

^<f^, (A3) 

which is the equivalent of equation If the large-scale magnetic field is weak enough so as to not significantly 
affect the dynamics in the main body of the disk (in particular, if thermal pressure supports the disk vertically 
against gravity), then pQ which appears in equation (IA3|) should be interpreted as the thermal pressure. However, the 
assumption of weak fields is not required, and equation (jlOp . which is derived from equation (jA3p . applies when the 
field is stron g as well as when it i s weak. Only in the parts of the main body of the paper where we combine equation 
^ with the iShakura fc SunvaevI (|1973D solution are we assuming that the large-scale magnetic field is dynamically 
weak. 

Equation (|A3[) is not the same as assuming that the magnetic and thermal energy densities are comparable at z « /i, 
which is sometimes quoted as the condition for suppression of the MRI. We note, however, that the condition on 
the Alfven speed is more fundamental, and the condition on the magnetic energy density is merely derived from it 
under a specific set of circumstances in the main body of the disk. In fact, analytical studies of both stratified and 
unstratified disks generally suggest that the condition k ■ va ^, is most important for MRI supp ression, regardless 
of the overall magnetic energy density (jBlaes fc Balbu£lll994l ; [Gammie fc BalbusllT99l iKim fc Os trikcr 2000). Thus, 
we believe that equation (jA3p is correct. 

In fact, what is most clear from the above MRI studies is that the stability criteria in each coordinate direction 
are roughly independent; for example, a strong toroidal field does not significantly affect the most unstable vertical 
wavelengths. Thus, if we were to suppose that the MRI is suppressed when the local magnetic and thermal energy 
densities are equal, it would be reasonable to assume that the appropriate condition is -B^^ > Sttp/j (i-e., that it involves 
the vertical rather than the total magnetic energy density), in which case we would derive (ph/po) ^5 Bzfi/SnpQ rather 
than equation (|A3|1 as the condition for MRI suppression; here, the subscripts have their usual meanings. We therefore 
see that the two possible conditions are ess entia lly the same for a disk that is dominated by gas pressure up to z « /i; 
one would simply need to modify equation (|A3[) by introducing a ratio of temperatures To/T^ on the right hand side, 
which is generally no more than a factor of a few. Only if the disk is radiation-dominated at z « /i is there a significant 
difference between the two conditions for suppression of MRI turbulence. In the case of a radiation-dominated disk, 
our use of equation (jA3[) means that we are assuming the MRI may be suppressed in the surface layer of the disk even 
when the magnetic pressure is weaker than the radiation pressure there. 
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B. MAGNETIC FIELD ADVECTION IN QUASI-STATIONARY ACCRETION DISKS 

Our work in this paper is concerned with accretion disks that are fundamentally time-dependent. In particular, 
we are studying situations in which a large-scale magnetic field is being dragged inward and causing the magnetic 
flux in the inner region of the disk to change in accordance with equation ([3]). Such disks are complicated to study 
analytically, and for this reason we used order of magnitude approximations to reach some of our conclusions in HA.'Il 
(in particular, to estimate Hb, the scale height of the vertical magnetic stress). 

In this section, we will give arguments for these approximations by considering the simplified case of a "quasi- 
stationary" disk. By this we mean a disk that is allowed to vary on the long (viscous) timescales at which magnetic 
field advection takes place but which is assumed to quickly adjust its structure on shorter timescales (in response to 
the changing magnetic flux) so that in most respects it can be treated as stationary. 

We begin by evaluating the (j) component of the induction equation ^ at z — h (where the magnetic diffusivity is 
negligible) in an axisymmetric disk, which gives 



dt 



[V4,Bz 

az 



v,B, 



V^Br) 



(Bl) 



If we assume locally circular orbits in Newtonian gravity and evaluate spatial derivatives of the orbital velocity to 
first order in /i/r, we can combine equation (jBip with the condition V • B = to derive 



dB, 



4>h 



dt 



'2^K ( Brh 



— iJzh 
r 



d d 



(B2) 



The first term in this equation represents the production of azimuthal field via Keplerian shear and is generally the 
most important effect; it leads to the creation of B^ on a characteristic timescale of ~ fl]^^ (i.e., the orbital timescale). 
Therefore, in accordance with our quasi-stationary approximation, we can assume that the disk quickly adjusts its 
value of Bcjih in response to the magnetic field advection so that dB^ph/dt is negligible on our timescales of interest. 
We can therefore rewrite equation (|B2() as 



dB, 



4>h 



dz 



'2~z 



K 



Br 



-B, 



B,, 



dlnv^h B. 



dz 



(t>h 

r 



Vrh 9 In (VrhB^h) 



Vzh 



dim 



(B3) 



We can simplify this equation by assuming that pv^ is constant with height so that there is no net buildup of mass in 
the vertical outflow. This is a reasonable assumption for thin disks for the timescales of interest. Thus, dlnvzh/dz — 
—din. ph/dz. Numerical simulations show that this quantity is, in turn, well-estimated by ~ H^^ throughout the 
atmosphere of the disk, even in time-d epende nt systems in whi ch radiation and magnetic fields help to provide vertical 
support against gravity (jHirose et al . 2006; K rolik et al.ll2007f) . Equation (|B3p then becomes 



dB. 



4,h 



dz 



Vzh 



B 



rh 



— Ozh 

r 



B. 



H 

H \ ^~ 



4>h 



Vrh 9 In (VrhB^h) 



Vzh 



dlnr 



(B4) 



If we combine this equation with the definition of Hb, use V • B = 0, and then ignore terms of order ^ H/r (assuming 
that the logarithmic radial derivatives are of order unity, which should be true except at boundary regions where the 
magnetic field structure changes dramatically), we obtain 



Hp 



H 



(B5) 



where Jb = {vzh/ HQ.k) {—B^h/ Brh) is typically a positive number in any magnetic field geometry (due to Keplerian 
shear). We therefore see that unless Vzh approaches the sound speed, we have Hb ^ H for typical generic field 
geometries. This corresponds to our approximation in the main body of the paper. Furthermore, since Vzh is simply 
the initial speed at which material is launched off the disk surface (before it enters any jet acceleration region), we do 
not expect it to be large; a reasonable estimate might be Vzh ~ {vr), where {vr) ~ ct {H/r) HVLk is the speed at which 
material in the main body of the disk accretes inward. In that case, we would have Hb H, which is even more 
favorable for inward advection. 

The above analysis shows that the vertical magnetic stress potentially available at z = /i in an accretion disk is 
more than enough to drag the surface layer inward. In fact, if we consider equation ()11|) in light of these results, the 
immediate suggestion is that the surface layer can advect inward much faster than the main body of the disk. The 
problem with this conclusion, however, is that the above analysis only considered the atmosphere on its own, without 
regard to the main body of the disk below it. A large shear between the main body of the disk and the surface 
layer is unlikely to be stable. Although the physics in the interface between the turbulent body of the disk and the 
nonturbulent region above it is complicated, we can estimate the effects of a large vertical shear in Vr in the case of 
ideal MHD. We use the r component of equation ([T]), which is 



dBrh d 
dt dz 



VzBr),^ , 



(B6) 



so that there is a term Bzhdvrh/dz which tends to produce Brh on a typical timescale (dvrh/dz) . Thus a large 
vertical shear in Vr will rapidly change the radial magnetic field — and thereby the vertical stress in accordance with 
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equation (|B4p — in the appropriate direction for the shear to be reduced. This suggests that very small values of Hb 
are not sustainable, and while it is difficult to predict the exact behavior, we expect that in the general case the 
base of the nonturbulent surface layer will advect inward at a similar speed as the main body of the disk. Note that 
turbulent stresses may also play a role in this region in ensuring that the disk and surface layer advect inward together 
(similar to a "friction" term). We do not consider their effect here except to note that they are probably smaller 
than the large-scale magnetic stresses. (Recall that the large-scale magnetic energy density is by definition com- 
parable to the thermal pressure at z = h, while the turbulent vertical stresses are likely to be smaller by a factor of ~ a.) 
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